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Abstract 

We determine the relative pion mass shift M V (L)/M n — 1 due to the finite spatial extent L 
of the box by means of two-flavor chiral perturbation theory and the one-particle Liischer 
formula. We use as input the expression for the infinite volume tttt forward scattering 
amplitude up to next-to-next-to-leading order and can therefore control the convergence of 
the chiral series. A comparison to the full leading order chiral expression for the pion mass 
in finite volume allows us to check the size of subleading terms in the large- L expansion. 



1 Introduction 



Shifts in particle masses and matrix elements due to the finite extent of the box are systematic 
effects in the Monte Carlo treatment of any lattice field theory. These shifts become particularly 
large when the spectrum contains light particles, and apply to all particles, no matter how heavy, 
provided they couple to the light ones. In QCD with light quarks, e.g. the nucleon mass receives 
such a correction, which disappears only as the product of the pion mass times the length of the 
box (M^L) gets large. 

Fortunately, as long as pions are light, chiral symmetry imposes strong constraints on the 
way observables deviate from their infinite-volume limit: chiral perturbation theory (CHPT) 
allows one to perform a systematic expansion around this limit and hence to control the finite- 
volume effects analytically [U 121 EJ El E] • in some cases one can directly translate the results of 
finite volume simulations into information about infinite- volume quantities without any need to 
extrapolate in the box size L jH] . CHPT is an expansion in the pion mass and particle momenta 
which have to be small in comparison to the chiral symmetry breaking scale, usually identified 
with 4nF n . The conditions of applicability read 

P «1, (1) 



4irF n ' 4ttF,. 

Particles inside a box of spatial length L with periodic boundary conditions and comparatively 
large or infinite time extent may only have discrete values of their spatial momenta, pk = 2Tmk/L 
with rik G Z. In this case the first condition in becomes a bound on the box size: 

£»^r~lfm. (2) 

Note that these conditions have no say whether L is large compared to the Compton wavelength 
of the pion or not. Both options are acceptable [3], but they imply different ways to organize 
the chiral series: 

M W L > 1 "p-expansion" (3) 

M W L < 1 «-> "e-expansion" . (4) 

We shall here restrict ourselves to the former case but investigate what the conditions (JTJd, |2j) 
mean quantitatively, by considering more than the leading order in the chiral expansion. The 
underlying assumption is that CHPT itself will tell us when the conditions of applicability are 
not respected any more, via a bad convergence behavior. 

In this paper we study the pion mass, M 7r (L), defined as eigenvalue of the QCD Hamiltonian 
in a LxLxL box (with periodic boundary conditions), as it is extracted on an Euclidean lattice 




Figure 1: The mass shift due to quantized momenta in the self-energy corrections amounts to a 
finite-size effect from pion exchange "around the world" (left), depicted by a "thermal insertion" 
(cross) in diagrammatic language (right). 
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for sufficiently large time T. The finite volume shift, M n (L) — M n , is - in coordinate space view - 
due to the possibility of the pion wrapping "around the world" or - in momentum space view - 
due to the discrete momenta, as depicted in Fig.^ The ultraviolet properties of the theory are 
untouched - with toroidal boundary conditions no new counterterms are introduced and the 
finite size effect is automatically finite. Our goal is to study the shift M n (L) — M n (with M % 
we always denote the infinite-volume value) as a function of M n and L. A preliminary account 
has been given in [7j. Here, we include one more order in the chiral expansion (at large L) and 
carefully analyze the uncertainty due to the errors of the QCD low-energy constants involved. 
In addition, we compare to the full one-loop result by Gasser and Leutwyler (which is leading 
in the chiral expansion, but includes the large- L suppressed terms). 



2 The pion mass in finite volume 
2.1 Finite volume calculations in CHPT 



In a series of remarkable papers [2l|4ll2], Gasser and Leutwyler have shown that chiral symmetry 
does severely constrain physical observables at low energy even if the system is confined to a 
finite box. The only condition is that the box length L must be several fermi, cf. (|2~|). 

Using the method of the effective chiral Lagrangian, one can establish low energy theorems 
for physical quantities of interest, which have the generic form 



Q 



Ml 



(5) 
(6) 



and the coefficients are quantities of order one in the chiral expansion. They depend on the 
low-energy constants (LEC), on ratios of momenta and quark masses, and, in finite volume, on 



A = M W L 



(7) 



For the pion mass and decay constant, the coefficient q\ in the p-expansion (JHJ) has been 
explicitly evaluated in Ref. 0, 



M n (L) 
F n (L) 



1 + 



2Ni 



-^i(A)+0(f) 



N, 



with 



= yz / dxe ~ 

Jo 



.|( n 2 +n 2 +n 2 )A 2 



(9) 



(10) 



where the sum runs over all integer values of n 12 ,3, excluding the term with (ni, n 2 , n 3 ) = (0, 0, 0) 
as indicated by the prime over the summation symbol. For a given value of n := n\ + n\ 
the integral can be performed analytically, and we can rewrite g\ as 



Tin 



n=l * 
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Table 1: The multiplicities m(n) in (jllj) for n < 20. 
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where K\ is a Bessel function of the second kind and the multiplicity m(n) indicates how many 
times the value n is generated in the sum in Eq. (jlUj) . The values of m(n) for n < 20 are given 
in Tabled For large argument the Bessel function K\ drops exponentially: Ki(x) ~ e~ x /*Jx. 
Since A ^> 1 by assumption, the sum (JTTJ) converges very rapidly, and it is easy to check how 
many terms are needed to get a good approximation for the complete sum. The convergence 
of the chiral expansion, on the other hand, is more difficult to test: the coefficient q2 has not 
yet been calculated neither for M n nor for F w . This would require a full two-loop calculation in 
CHPT in finite volume. We are going to argue, however, that there is a fast and reliable way to 
check the convergence of the chiral expansion - a way that is based on a formula due to Liischer 
|H| which we will discuss in detail. Before doing so, it is useful to rearrange Eq. (JHHHJ) by making 
explicit the expansion of each coefficient q m in a series of rapidly decreasing exponentials (in A) 



il4(L) 



Mir 



Mir 



i + (egiY+e 2 ^+o(e 3 ))^(A) 

M 



F n (L) 



+ {^12 + + 0(e)) MV2X) + ... 

oo oo 

i+EE£ m ^i(V^A) 

m=l n=l 

oo oo 



(12) 
(13) 



m=l n=l 

where at leading order in the chiral expansion the coefficients q± n follow from Table 1, 

2m(n) v 2Nfm(n) 



Hln 



Qln 



n A 



(14) 



Notice that the Bessel functions that we have factored out in Eqs. ()12[ IT3*j) need not necessarily 
appear in this form also at higher orders in the chiral expansion - the coefficients q^^f are 
presumably complicated functions of A, in general. We do expect, however, that also at higher 
orders the result may be expressed as a series of exponentials, whose leading asymptotic behavior 
should be captured by the Bessel functions in (jl2|13|) . 



2.2 Liischer's formula for M 7r (L) 

Liischer considered the problem of evaluating the finite size effects on a particle mass from a 
different point of view. Via graph-theoretic arguments he proved an elegant relation between the 
(Euclidean) finite volume mass shift and the (Minkowski space) tttt forward scattering amplitude 
F(v) in infinite volume [Bj: 

o /*oo 

i\4(L) -M 7T = -J^r^-jr J dyF(iy)e~V^t + (e- ML ) . (15) 

The integration runs along the imaginary axis, where the amplitude is far away from its cuts (cf. 
Fig.EJ): only the real part of F(iy) contributes to the integral. The one-particle Liischer formula 
(I15|) is an asymptotic relation for large L; it is proven to all orders in perturbation theory for 
an arbitrary massive QFT, and the subleading piece is tamed by the bound M > ^/3/2 M n 
jS]. An extra term originating from the 3-particle vertex in the original formula jH] does not 
appear in our case, due to the odd intrinsic parity of the pion and parity conservation in QCD. 
Since eq. ()15j) builds on the unitarity of the theory, it holds only in full QCD. A variant for the 
quenched approximation is not known, if at all possible. 

Neither the formula, nor its derivation make use of the chiral expansion (indeed the for- 
mula was derived before CHPT was applied to finite volume effects). However, for practical 
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Figure 2: Integration contour in the complex v plane where v is the crossing variable in the 
(Minkowski space) forward scattering amplitude. The cut is generated only at NLO in the chiral 
expansion. 



applications we need to insert an explicit, analytic representation for the tttt forward scattering 
amplitude, and for this we shall rely on CHPT. A comparison to the expression (jl2j) shows that 
by inserting the chiral expansion for F in and evaluating the integral one can obtain the 
coefficients q^ x . If we write (extracting a factor 16-7T 2 for later convenience) 

F(iy) = 16tt 2 [£F 2 + e 2 F 4 + £ 3 F 6 + 0(£ 4 )] , (16) 



the integration of the term F 2m yields the coefficient q^ x . Since for 2 degenerate dynamical 
flavors the tttt scattering amplitude is known to next-to-next-to-leading order [5| in the chiral 
expansion, i.e. up to F 6 , we can evaluate the coefficients q^ 1: m = 1, 2, 3 for Nf = 2. All this 
without having to do a single new loop calculation in CHPT in finite volume. 

We stress that the Luscher formula is designed to yield the leading term in the large L 
expansion: whenever the relative suppression factor 

subleading = / _ { M_ Mit)L) \ ,^ 
leading V / 

is not small, subleading effects may be of practical relevance. In order to estimate subleading 
effects in the large L expansion we will go back to the expansion of (?i(A). 

The rest of this article is devoted to such a combined use of Luscher's formula (|15jl and the 
chiral expansion (jlfijl for two-flavor QCD. It is useful to keep in mind that the chiral formulae 
contain low-energy constants that have been pinned down in real Minkowski space experiments; 
some stem indeed from tttt scattering data, which determine the amplitude in the physical region, 
e.g. for forward kinematics, on the cuts in Fig.|21 The very fact that the chiral expression for 
F(y) holds in the whole complex v plane is a key ingredient of this work - we rely on the good 
analyticity properties of the chiral amplitudes. 

2.3 Combining CHPT and the Luscher formula 

The forward scattering amplitude which is needed in Luscher's formula can be expressed in terms 
of the isospin invariant amplitude A(s, t, u) as follows 

F(u) = T(2M 7r (M w + u),0,2M w (M n -u)) 
T(s,t,u) = A(s,t,u) +3A(t,s,u) + A(u,s,t) , (18) 

where the only independent kinematical variable v reads 

v = = -i_ - M w , (19) 

4 



in terms of the momenta p a ^ of the two initial pions. Since A(s,t,u) is symmetric under t, u 
exchange, the forward scattering amplitude is an even function of v: F(i/)—F(—u). 
We find it convenient to use a dimensionless integration variable and introduce 

y = y/M 7T , v = v/M w . (20) 

In these variables and with N=16tc 2 the chiral expansion of the forward scattering amplitude is 

F{u) = N [iF 2 {u) + £ 2 F 4 (i>) + £ 3 F 6 (i>) + 0(£ 4 )] . (21) 

If we combine Liischer's formula with the chiral expansion for F we obtain a simple expression: 

M:M.)-M:: 



R M (M n ,L) 



= R M (M n ,L) 

3 ro ° 
A 



/OO 
dy e~^ x [£F 2 (iy) + £ 2 F 4 (iy) + £ 3 F 6 (iy) + 0(£ 4 )] 
-OO 



+ 



j [e/ 2 (A) + eh(\) + e 3 / 6 (A) + o(e)] + . . . 



(22) 



where the ellipsis indicate terms of order e A/,rL and 

/OO 
dye-^ x F 2m {iy) . 
■OO 



(23) 



Using the expression for A in [H] (which is specific to Nf — 2) and splitting its 6j coefficients 

= 6;/iV for i = l,...,4 
= 6;/iV 2 _ for i = 5,6 
= 6? + eR for z = l,. ..,4 , 



(24) 



the calculation of the coefficients l^m in 
7 2 (A) = -S°(A) 



can be made in large parts analytically: 
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(25) 



where the integrals B 2k admit a simple analytical representation 



/CO 
dy f k e -^ x 
-oo 



r(fc + i/2) (r 
r(3/2) [x) Kk+l{x) 



B\\) 
B 2 (X) 



2 MX) , 

\k 2 {x) . 



(26) 



The integrals denoted by R k , i — 0, . . . , 4 are defined as follows 

J2?(A) = iV 2 |^ jH rfy fe'^ x Kr (2(1 + ijr)) for j 



even 
odd 



(27) 



(the functions appear in the im scattering amplitude at the two-loop level P - there 
they are defined without the nir superscript - and collected in Appendix B) and can be easily 
calculated numerically, while an analytic representation is not available. The integrals B 2k 
and R k are plotted in Fig-ISl As can be seen there, they are all of a similar magnitude - in 
the representation (}2*2*|) the small parameter £ appears explicitly such that all the remaining 
coefficients are not expected to show any special hierarchy. 

The last ingredient needed for an evaluation of Rm(M^, L) is the numerical values of the 
LEC which appear in the coefficients b{ (which are specific for Nf — 2), and in F % , if the latter 
needs to be computed from M n . The expression of the hi in terms of the relevant LEC can be 
found in the appendix. At order p 2 no LEC appears. At order p 4 four LEC appear: three of 
them were determined rather precisely in JU], whereas for the one that dictates the NLO quark 
mass dependence of M 2 , we rely on the estimate in pQ. Altogether, this means that we use 



Ax 
A 3 



fi 2 

U.l/ -0.03 



0.59 



1-1.40 
-0.41 



GeV 
GeV 



A 2 = 1.20 +° f 6 GeV 
A 4 = 1.25 ±8:il GeV 



(2f 



At order p 6 six new LEC appear. Two of them were determined in ^Hj, whereas for the 
remaining four we rely on the resonance saturation hypothesis with a (conservative) 100% 
error estimate: 



f x = -1.5x(l±l) 
f 3 = -4.2 x (1 ± 1) 
f 5 = 3.8 ±1.0 , 



f 2 = 3.2 x (1 ± 1) , 
f 4 = -2.5 x (1 ± 1) 
f 6 = 1.0 ±0.1 . 



(29) 



Inserting these values and evaluating the corresponding bi one can then determine the integrals 
I 2 , Iq in the ^-expansion (J22J), and hence RM(M n , L), the details being given in section 3 and 
the appendix. 



3 Numerical evaluation 

3.1 Quark mass dependence of F n 

In CHPT the expansion parameter is £, and in (fT2*l |2*2*|) the finite volume effects have been 
expressed as a power series in this parameter. If, in a lattice calculation, both M 7r (L) and 
F,r(L) have been determined, then the square of M 7r (L)/(47rF vr (L)) may be taken as a first 
approximation to £, with iterative refinement through (j!2[ [T3|) . In order to give a numerical 
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Figure 3: The integrals B 2k and R k which are needed for the evaluation of the finite volume 
corrections to NNLO. The B 2k are known analytically, the R k have been determined numerically. 



prediction with only M n as input, we must know how £ depends on the quark mass, or - 
equivalently - we are left with the problem of evaluating the pion mass dependence of F^. The 
analytic expression of the latter is known to next-to-next-to leading order ^T], and reads in 
our variables 



Fll + X 



L + L 



--L 2 + L 
4 



6 



3 



29 \ 
12/ 



13 
192 



(30) 



where L = log(/i 2 /M^), X = M 2 /(NF 2 ) and f F (/j) is the relevant combination of 0(p 6 ) LEC. 
In order get a numerical value of F^ at fixed we need to specify our input parameters; the 
£i have been given in the previous section, whereas for the new LEC we take fp(fi) = ± 3 and 
vary the renormalization scale between fi = 0.5 and 1 GeV. The only remaining parameter is 
the value of F, the decay constant in the chiral limit. We fix it by inverting Eq. ()30|) - now at 
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Figure 4: The plot on the left shows the lex-band of the pion mass dependence of F„. The plot 
on the right shows the dependence of £ on the pion mass, without including any uncertainty. 



the physical pion mass - and expressing F in terms of F n . From F„ = 92.4 ± 0.3 MeV we obtain 

F= (86.2±0.5)MeV . (31) 

The pion mass dependence of F n is illustrated in Fig.0] and is seen to be rather mild. This 
means that in the finite-volume condition (J2J) the numerical value on the r.h.s. (which was 
assigned using the physical pion mass) holds to a good approximation also for substantially 
heavier pions - the increase of F n with the pion mass is too mild to matter in this respect. 
In our numerical studies of the CHPT formulae we restrict ourselves to L> 1.5 fm and actually 
expect that this bound might already be too low with respect to the condition (J2J). Let us remark 
that the same bound has to be respected even if one works in the e-regime (j3J and wants to 
compare to CHPT formulae: in this regime one must still have a relatively large volume and is 
consequently forced to use tiny - smaller than physical - quark masses. 

By contrast, the parameter £ remains small even for pion masses of about half a GeV, as 
is shown on the right-hand side of Fig.HJ In our numerical analysis we will use £ exactly as 
given in this figure and ignore the uncertainties of this "computed" £ (shown only in the plot on 
the left-hand side of Fig.HJ) since, as explained above, by measuring F n (L) one can determine £ 
iteratively from the lattice data. We do, however, determine the uncertainty of the coefficients 
I<im in the series (|22|) , and for that aim we use the correlation matrix among the input parameters 
which can be obtained from Ref. JU|- The numerical values and the details of the analysis are 
given in Appendix C. 



3.2 Finite volume effects in the pion mass 

We are now in a position to evaluate the Liischer formula (}2"2"j) for the relative mass shift 
Rm^M^, L). Before doing so, we find it instructive to have a look at the product 

F{iy)e^^+y L (32) 

which is the integrand in the Liischer formula 1)15)1 . 

We have plotted this integrand in Fig. El for different values of M„ and L, evaluated at LO, 
NLO and NNLO in the chiral expansion. Since the calculation is based on CHPT and the 
integration variable y has the dimension of a mass, the integrand must tend to zero sufficiently 
fast beyond O(100)MeV - otherwise the outcome of the calculation cannot be trusted. This 
criterion results in a clear veto against too small box sizes (around ~ lfm). By the absolute 
amount, the two integrands are close to each other whenever the pion is light and the box is 
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Figure 5: The integrand I(y) = e~V^+^ L F(iy) with LO, NLO, NNLO input from CHPT. 
Note the change in scale in the lower panel. 



sufficiently large. However, in the whole range of and L shown in the graphs, the relative 
difference between the LO and the NLO integrand is substantial. The reason is that at LO the 
function F(u) is a constant (F 2 can depend at most linearly on u, and since it has to be even 
in v it can only be a constant): the energy dependence of the forward scattering amplitude - in 
particular the one dictated by the unitarity cuts - appears only at NLO. The difference between 
LO and NLO is therefore not very useful in judging the convergence of the chiral series, which 
should rather be evaluated by comparing NLO and NNLO. The plots in Fig. 03 show that even 
for M n = 500 MeV the correction in going from NLO to NNLO is not unacceptably large. Note 
that in the graphs in Fig.|3]the product M n L ranges from about 1 (top left) to 10 (bottom right): 
in the former case the product is too small for the asymptotic formula (JT5j) to apply - we will 
confirm this below, by comparing this result to the one obtained with one-loop CHPT in finite 
volume, eq. (jHJ. 

The relative mass shift Rm{M^^L) calculated with the Luscher formula ()15 [ 122 (1 is plotted 
in Fig. El as a function of L for different values of M n and in Fig. [7| as a function of M n for 
selected values of L. In the latter figure we also show the uncertainty coming from the LEC - 
as expected the band grows with the pion mass. These figures show that the chiral expansion 
converges nicely for very light pions (M n = 100 MeV) and more slowly, but still satisfactorily 
for heavier pions. For light pions, however, the comparison between the full one-loop CHPT 
calculation and Liischer's formula (evaluated with the LO forward scattering amplitude) shows 
that the leading exponential behavior is not numerically dominating, even for volumes as large as 
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Figure 6: The relative mass shift Rm{M v , L) versus L for a few values of the (infinite volume) 
mass M n . In addition, the full 1-loop result (JBJ) by Gasser and Leutwyler is included. 

L = 4 fm (where M n L~2 for the lightest pion considered). Here, the use of Liischer's formula is 
not justified, and one should rather rely on the full one-loop result in CHPT in finite volume. On 
the other hand, the fact that the NLO correction in Liischer's formula is rather large indicates 
that even the full one-loop result does not give a reliable answer. In such cases (e.g. when one will 
be able to simulate 200 MeV pions in a 2 fm volume) one would need a full two-loop calculation 
of the pion mass in CHPT in finite volume in order to reliably estimate the finite size corrections. 

For heavier pion masses the higher exponentials which are neglected in Liischer's formula 
are less important and one is entitled to fully rely on this convenient formula. We find that 
for masses above 200 MeV and L > 2 fm the finite size effects are at most of the order of a 
few percent. We stress that in this range of masses and volumes our evaluation of the finite 
size effects is reliable because we are able to check both the convergence of the chiral expansion 
and that of the large volume expansion of Luscher. For M n L ~ 4 and L ~ 4 fm or larger the 
predicted shift is below 0.1%: in such cases we may conclude that for all practical purposes the 
mass obtained on the lattice coincides numerically with the infinite-volume one. Mass shifts of 
the order of 1% are predicted only for M n ~ 200 -j- 300 MeV, L ~ 2 2.5 fm). This is the region 
where precision tests of lattice QCD calculations (we stress again that our calculation applies 
only to full QCD) will need the application of such corrections. 

Our numerical findings are summarized in Tab. El where we give for selected values of 
and L the relative mass shift Rm computed via the Luscher formula (fTK| l2"2"j) with LO, NLO 
and NNLO input from CHPT. For comparison, we give the full one-loop result without large L 
expansion (fT2l H3J) , due to Gasser and Leutwyler [2] . We have also combined the two results by 
adding to the NNLO Luscher formula the series of large- L suppressed exponentials which appears 
in the full one-loop result but not in Liischer's formula. This last figure is our best estimate 
of the total finite-size correction, but a sizable difference between this and the NNLO Luscher 
result signals the presence of large uncertainties. In such cases (as we have argued above) a full 
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Figure 7: The relative mass shift R^(M n , L) versus M w for a few values of the box length L. The 
thin lines around the NNLO curves represent the corresponding uncertainties. For comparison, 
the full one-loop result (jHJ) of Gasser and Leutwyler is included. 

two-loop evaluation of these corrections would presumably settle their value. 

Note that in the first column of Tab. El we give the finite size effects for L — 1.5 fm, although 
for such small volumes we are very likely already in conflict with @, and are therefore outside the 
region of applicability of CHPT. We include this column in the hope that future high-precision 
lattice data in the regime 1.5fm<L<2fm might pinpoint what eqn. (J2J) means quantitatively. 

4 Summary and Discussion 

Making all the necessary extrapolations for lattice QCD calculations (lattice spacing to zero, 
quark masses to their physical value and volume to infinity) at the same time would be enor- 
mously expensive in terms of computer time and is practically unfeasible. Wherever possible one 
should try to use analytical methods as an aid. As far as the finite-volume effects are concerned, 
the necessary theoretical tools to control these artifacts have been developed in the eighties by 
Liischer |Hj and by Gasser and Leutwyler 011]. 

Recent progress on the side of lattice calculations make now this issue of concrete relevance, 
and gave us the motivation to make a thorough numerical study of the finite size effects on the 
pion mass. In order to do so we have explicitly evaluated the Liischer formula using as input the 
forward scattering amplitude evaluated at leading, next-to-leading and next-to-next-to-leading 
order in two-flavor CHPT. Pushing the chiral expansion to such a high order has shown to be 
absolutely essential in order to obtain good control on the convergence of the chiral series: for a 
wide range of pion masses and lattice volumes we find that the we can evaluate these corrections 
reliably. As shown first by Liischer, these corrections vanish exponentially with M n L and are 
therefore negligible for sufficiently large masses and/or volumes. We have determined the regions 
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3.31 


2.73 x 10~ 2 

5.15(14) 

5.60(14) 

8.74 

11.6 


1.46 x 10~ 2 

2.41(7) 

2.58(8) 

3.85 

4.97 


0.83 x 10~ 2 

1.26(4) 

1.34(5) 

1.91 

2.41 


5.02 x 10~ 3 

7.15(24) 

7.51(29) 

10.2 

12.7 


3.13 x 10~ 3 

4.27(15) 

4.45(18) 

5.81 

7.13 


200 MeV 


2.87 x 10" 2 

6.22(47) 

7.61(73) 

6.56 

11.3 


1.08 x 10~ 2 

2.03(18) 

2.39(27) 

2.00 

3.31 


0.45 x 10~ 2 

0.78(8) 

0.90(12) 

0.73 

1.17 


2.04 x 10~ 3 

3.31(34) 

3.74(52) 

2.95 

4.65 


0.96 x 10~ 3 

1.50(16) 

1.67(24) 

1.29 

2.00 


0.47 x 10~ 3 

0.71(8) 

0.78(12) 

0.59 

0.91 


300 MeV 


1.31 x 10~ 2 

2.87(41) 

3.65(77) 

2.25 

4.59 


0.39 x 10~ 2 

0.78(12) 

0.95(22) 

0.56 

1.12 


1.27 x 10~ 3 

2.42(41) 

2.89(72) 

1.65 

3.27 


0.45 x 10~ 3 

0.82(14) 

0.96(25) 

0.54 

1.05 


1.64 x 10" 4 

2.93(53) 

3.41(92) 

1.88 

3.65 


0.62 x 10" 4 

1.09(20) 

1.26(35) 

0.69 

1.32 


400 MeV 


0.58 x 10" 2 

1.35(28) 

1.75(55) 

0.85 

2.01 


0.13 x 10~ 2 

0.30(6) 

0.37(12) 

0.17 

0.41 


0.34 x 10~ 3 

0.74(17) 

0.91(32) 

0.40 

0.97 


0.94 x 10" 4 

1.98(46) 

2.41(86) 

1.03 

2.50 


0.27 x 10" 4 

0.56(13) 

0.67(25) 

0.29 

0.69 


0.79 x 10~ 5 

1.64(39) 

1.96(73) 

0.83 

1.99 


500 MeV 


0.26 x 10" 2 

0.66(17) 

0.87(34) 

0.34 

0.95 


0.47 x 10" 3 

1.16(31) 

1.50(60) 

0.54 

1.58 


0.09 x 10~ 3 

0.23(6) 

0.29(12) 

0.10 

0.30 


0.20 x 10" 4 

0.48(13) 

0.61(26) 

0.21 

0.62 


0.44 x 10" 5 

1.07(29) 

1.34(57) 

0.45 

1.35 


0.10 x 10~ 5 

0.24(7) 

0.31(13) 

0.10 

0.31 



Table 2: The relative finite size effect Ri^iM^, L) for selected values of and L. In each cell we 
give the result of the Liischer formula (JT5j) with the forward amplitude at LO, NLO and NNLO 
accuracy from CHPT. The fourth entry gives the GL result (jHl fTTj) and the last one contains the 
full one-loop result shifted by the difference between the NNLO and the LO Liischer formula - 
this, we believe, is a reasonable estimate of the total effect. The power of ten given in the first 
entry of the cell applies to all other entries. Note that the first column is very likely in conflict 
with the condition (J2J). Note, finally, that the GL result fTTj) for the relative shift of F w is 
(—4) times the fourth entry (for Nf =2). 

where these corrections are important for precision calculations. Our numerical results are given 
in Figs. El [7|and in Table |2] these are the main results of this paper. 

If M n L is not very large, keeping only the leading exponential in the large volume expansion, 
as is the case in Luscher's formula (|15j) . may not be accurate enough. In such cases one can 
take into account the whole series of exponentials by working with CHPT in finite volume. One- 
loop expressions for the pion mass and decay constant are available in the literature [2j and, 
in numerical form, in our Table El and we have used them to estimate numerically for which 
values of M n L is Luscher's formula not enough accurate. For such situations, it is at the moment 
difficult to estimate the size of finite- volume effects, and a full two-loop calculation in CHPT in 
finite volume seems the only viable way to do it reliably. 
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A: The coefficients 61, ... , b$ 

For convenience we provide in this appendix the explicit expressions for the effective coupling 
constants 4, . . . ,6 6 which appear in the Tin scattering amplitude up to order p 6 jHJ HD] in the 
split convention 



where 



7 ~ 4 ~ 1 ~ ~ 13 
1 6 3 2 A 18 

b\ = --L 2 + {--h- —i 2 - 4 - —h -— U + -44--4 2 -3£ 3 i 4 -5£i 

1 g I 9 9 3 108J 3 2 3 3 4 

28- 80- 15- 26- 34tt 2 3509 _ 

+ 27 fl + 27 £2 -T 4 + ¥^ — + 1296 +ri 
2 - 4 - 2 
bl = -L--4+24-g 

2 36 ^ l 9 2 2 3 3 54 I 3 4 

- 166- 9- 8- 317tt 2 1789 _ 

-44 4 + -4 4 H h r 2 

27 2 3 9 216 432 



1 - 1 - 1 - 7 
3 2 3 6 12 

11 40 ?2 f 38 « 20 « „; 365^1 ~ 4~~ 2~~ 

6; = L + <^ 4 £ 2 + 2 4 H \l + -44 + -44 

3 9 1 9 3 216/ 3 3 

89- 38- 7- 311vr 2 7063 „ 

H 4 H 4 4 1 V r 3 

27 9 3 432 864 

Y 85- 2 r7~ 107- 625 ^1 ~ 31- 145- 7vr 2 66029 

4 = — L + < -4 H 4 ^ fx 4 H h r 5 

72 1 8 72 288 J 36 108 54 20736 

T 5- 2 r5~ 25- 257^1- 7- 35- vr 2 11375 

4 = — L + < —4 H 4 >L 4 4 H h r 6 33 

24 1 72 72 864/ 108 108 27 20736 



-4 47 U+ 2 4 

! 216/ 3 


17tt 2 
216 


1655 

i h r 4 

2592 


625i 


l-i - 511, - 

r 36 


' 288 J 


257 1 


^ 7 4 


~ 864 J 


1 108 



A 2 2 

= log -f- , r t = N 2 rl(n) , L = log . (34) 



Note that the quark masses exclusively enter through £ and L; the remaining quantities are 
independent thereof. 
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B: T_he functions K 



7T7T 



The functions Kf % (x), i = 0, . . . , 4 have been introduced in jH], and we reproduce them here for 
convenience. With 

and h(x) = Ar j=. log ■ 



they read 



4 

1 - - 

x 



K™{x) 
K^(x) 



Kr{x) 



KF(x) 



Z 1 / N 2 

N k{x) + W 
zh 2 (x) 

ft.,, > 7T 2 , , . 7T 2 

+ -/»(*) - ^ 

1 

xz 



/CM 



6)x 



C: Evaluation of the uncertainties 



In order to quantify the uncertainty of the l 2m in (j22|- we need to know the correlation matrix 
among the LEC involved and the partial derivatives of the l 2m w.r.t the LEC. 

The correlation matrix CV, among our ten input parameters is given in Tab. El It has been 
obtained from Ref. JU|, which represents so far the best determination of the LEC appearing in 
the 7T7T scattering amplitude. We remark that some of our input parameters were also used as 
input in ^0] and therefore are statistically independent - this is seen in Tab. El where some of 
the off-diagonal matrix elements are zero. 



Combining and one finds 



dh 

W 1 



4B C 



3 ' 



dh = 8 o _ 32 2 dh 
' J 3 



d£ 2 



-B° 
2 



dh 

W 4 



-2£? c 



(35) 



and 



dh 

all 

dh 

d! 2 



3 3 4 9 



27 



560 ~ 128 ~ 

1j — 



80 







80 



128 



16 



-tin — tin 



4- & 



dh 

dh 
dh 
dl 4 

dh 
Of i 
dh 
dr 2 
dh 
dr 3 



15L - 5£ 3 - 1U 4 - - ) B u - 

4 / 



114 



104 

448 



5 U 



160 ~ 32 ~ 128 ~ 

L 1\ £ 2 + 

3 3 3 9 



110 

D 2 52 D o 64 D i 



55°(A) 
45°(A) 

85°(A) - 8B 2 (X) 



160 
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^ = 8£? (A)-565 2 (A) 
ar 4 

^ = 16B°(A)-48£ 2 (A) 
dr 5 

^ = 165°(A) + 16B 2 (X) . (36) 
or 6 

These two vectors are then used to sandwich the correlation matrix of our input parameters. 
With x = {£i, . . . , £±, fi, . . . , f 6 }, the resulting uncertainty is 



AR 



M 



(37) 



where dl^/dxi is zero for z > 4, and the rationale for omitting the uncertainty in our evaluation 
of £ has been explained in the text. 
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T2 



r 3 



r 4 



r 5 



r 6 



h 

h 
h 

h 



3.5 10- 1 -3.3 10- 2 -3.0 10- 2 6.7 1(T 2 -8.8 1.7 1(T 2 -1.2 -2.1 10" 1 -5.4 10" 1 -3.7 10' 2 

1.2 1(T 2 -3.0 10" 3 -7.2 10~ 3 7.2 10~ 6 3.7 10~ 3 9.8 10~ 2 -3.9 10" 1 1.1 NT 2 -4.6 10~ 3 



5.8 -5.5 10" 2 



-3.3 10" 2 7.0 10" 4 



4.8 10~ 2 -8.6 10" 6 2.3 10~ 3 -1.2 10 1 -7.4 10~ 2 -9.1 10 2 -2.210" 



2.3 



10 



18 



4.8 10" 4 -5.7 10^ 5 



- -2.1 1Q- 2 -5.8 10" 4 



2.3 1.6 HT 1 



6.3 



1.9 3.7 10- 1 



1.1 9.2 10" 2 



1.1 ID -2 



Table 3: The correlation matrix of our 10 input parameters. 



